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by 
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Summary 


A rotor ' section aerodynamic design package is developed by coupling Constrained 
Direct Iterative Surface Curvature (CDISC) design method with FPX rotor code. The 
coupling between the CDISC design and the FPX flow analysis is fully automated. The 
CDISC design method employs a predictor-corrector procedure iteratively to determine 
a surface geometry which produces a target pressure distribution, where target pressure 
distribution is either pre-defined or automatically generated through flow and geometry 
constraints. The FPX code is an extended Full- Potential rotor Computational Fluid 
Dynamics (CFD) code, which solves three-dimensional unsteady full-potential equation 
in a strong conservative form using an implicit approximate-factorization finite-difference 
scheme with entropy and viscosity corrections. Application of the CDISC design method 
coupled with the FPX rotor code is made for rotor blades in hovering motions. Several 
design examples are presented to demonstrate the capability and efficiency of the new 
package in rotor section design. 
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1. Introduction 


The recent advances in Computational Fluid Dynamics (CFD) have provided accurate 
and detailed flow solutions in an efficient way. These advances have generated interest in 
integrating CFD code into design methods. Automated aerodynamic design methods are 
beginning to become part of this process. These methods can be divided into two general 
categories: inverse methods that use an inverse solver, and predictor-corrector methods 
that use a direct flow analysis in an iterative manner. The predictor- corrector methods can 
be coupled with any flow solver and hence take advantage of the advances in Computational 
Fluid Dynamics. 

The Constrained Direct Iterative Surface Curvature (CDISC) design method of Camp- 
bell and Smith 1 and Campbell 2 is a pressure based predictor-corrector design method. The 
CDISC method determines a surface geometry which produces a target pressure distribu- 
tion, where taxget pressure distribution is either pre-defined or automatically generated 
through flow and geometry constraints. The method has been coupled with several ad- 
vanced CFD flow solvers 1-5 for fixed-wing section design, including design in the presence 
of nacelles. Results from these codes indicate that the CDISC method is robust and 
accurate. 

This report presents the coupling of the CDISC method with a rotor Computational 
Fluid Dynamics code - the FPX code. As the best of author’s knowledge, this is the 
first time on using the CDISC method in helicopter rotor section design. In the next few 
sections, a brief description of the CDISC method and the FPX rotor CFD code is given, 
which is followed by the implementation of their coupling. Finally, several design examples 
axe presented to demonstrate the capability and efficiency of the new design package in 
rotor section design. 

2. The CDISC Method 

The Constrained DISC method (CDISC method) is an extension of the basic DISC 
method. The DISC method employs a pressure based predictor-corrector design procedure, 
where the calculated surface pressure distribution is compared with a target distribution 
and the correction to surface geometry is then made to achieve the target pressure dis- 
tribution. In more detail, the DISC method starts with an initial surface geometry and 
flow condition at design point, a flow analysis code is used to generate surface pressure 
distribution. The calculated pressure distributions are then compared with the target dis- 
tribution and the difference gives the correction to the surface geometry through a design 
module. A new geometry is obtained and the flow analysis code is used again to generate 
a new surface pressure distribution until the target pressure distribution is achieved. 


The DISC design module relates the difference between the calculated and the target 
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surface pressure coefficients (A C p ) to the change in surface curvature (AC) at a given 
chordwise location. Different modules are used in different regions during the design pro- 
cess. For subsonic and low supercritical flows where a small change in surface geometry 
can only result in a small change in surface pressure, the relationship is derived 1 from the 
momentum equation as 

AC = A{\ + C 2 ) B AC V (1) 

where A is 1 for the upper surface and -1 for the lower surface; B is input constant 
ranging from 0.0 to 0.5, with the higher values giving faster convergence but decreasing 
the numerical stability; and C is streamline curvature at aerodynamic surface. For the 
moderate to high supercritical flows where a small change in surface geometry can cause a 
large change in surface pressure coefficients, the relationship is obtained 1 from supersonic 
thin airfoil equation as 

, „ d(AC r ) 1 (2) 

dx 2 [1 + (dy/dx) 2 } 1 5 


where y defines airfoil surface. 

The advantage of the DISC method over the inverse method is that the method can 
be treated as a “black box”. The “black box” can be incorporated into any CFD code, 
and thus take the advantage of the advances of CFD. 

The Constrained DISC method specifies desirable characteristics of the target pressure 
distribution. The target pressure is automatically generated and updated at each step 
based on the specified characteristics - the constraints. The typical flow constraints are lift, 
drag, pitching moment and pressure gradient. The typical geometry constraints are airfoil 
maximum/minimum thickness and leading edge radius, and so on. Including geometry 
constraints in the design process will yield an airfoil/ wing/ rotor blade that not only satisfies 
target pressure distributions but also be practical to build. In the CDISC design method, 
constraints are easily included in the design process to satisfy various design requirements. 
For example, upper and lower surfaces can be designed individually and therefore designing 
only one surface is possible, which is very useful in design of very thin airfoil sections. 
Partial surface design is an option. Surface smoothing eliminates sharp changes in surface 
curvature. The limitation of surface curvature is useful in highly cambered airfoil section 
design. 

The three-dimensional (3D) wing/rotor blade design is considered as a set of airfoil 
section design problems at specified span-locations. For 3D span-load distribution design, 
the DISC method can also be used by modifying twist distributions. For simultaneous 
airfoil section and span-load distribution design, the surface curvature and twist modifi- 
cation methods can be used together. The CDISC design method has been shown to be 
comprehensive and easy to use, it is therefore highly desirable to use the method in various 
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configuration design problems. In the present work, the CDISC is used for helicopter rotor 
section design when coupled with a rotor CFD code. 


3. The FPX Rotor Code 


While in the fixed-wing aerodynamic computational community more and more ex- 
pensive and complex Euler and Navier-Stokes methods are used recently, potential methods 
still serve as a major analysis tool in the rotary-wing aerodynamic computational commu- 
nity. The eXtended Full-Potential (FPX) rotor code 6 is o ne such accurate and efficient 
potential method, which represents an industry standard for rotary-wing computations. 
The FPX code is a modified and enhanced version of Full-Potential Rotor (FPR) code 7 . 
The code (either FPR or FPX) has been used in various helicopter hover and forward flight 
cases, including Blade- Vortex Interaction (BVI) calculations 8 , nonlinear acoustic analysis 9 
and coupling with the comprehensive helicopter code CAMRAD 10 . The application of the 
code produces excellent results. The code is also highly optimized and a typical steady 
run takes 1-7 minutes on a Cray-C90 supercomputer with a single processor. 


The FPX/FPR codes solve the unsteady three-dimensional full-potential equation in 
a strong conservation form using an implicit finite-difference scheme for flow around rotor 
blade. It is less expensive than either Euler or Navier-Stokes methods, and yet produces 
accurate solutions for various helicopter flows without significant separations. 

The unsteady, three-dimensional full-potential equation in strong conservation form 
in blade-fixed body-conforming coordinates (£, 77 , £,t) is written as 


d 

dr 



1„(pJL\ + 1.(PL\ + !L(PK 

dC J ' + drj K J ’ dC J 


) = 0 


(3) 


with -- .... 

p = {1 + ^-h-24, - (V + {<)$< - ( V + rn)i, ~(W + c,)*(] } * (4) 

where $ is the velocity potential, U, V and W are contravariant velocity components, p 
is the density, and J is the grid Jacobian. 

The FPX/FPR codes solve Eq. (3) using an implicit finite-difference scheme, where 
the time-derivative is replaced by a first-order backward differencing and the spatial- 
derivatives are replaced by second-order central differencing. The resulting difference 
equation is approximately factored into three operators X TJ , and in £, 77 and £ 
directions, respectively 

I € I,I c ($ n+1 - $ n ) = RHS (5) 

The detail of the scheme is presented in References 6 and 7. 
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The FPX is the substantially modified version of the FPR code. Both entropy and 
viscosity corrections are included in the FPX code. The entropy correction potential formu- 
lation accounts for the shock produced entropy to enhance physical modeling capabilities 
for strong shock cases. Either a two-dimensional or a three-dimensional boundary layer 
model is coupled with the FPX code to account for viscosity effects. In addition, an axial 
flow capability is added into the FPX code to treat tilt-rotors in forward flight. In addi- 
tion to the 0-H grid topology, an H-H grid topology is added as well. More recently, the 
Vorticity Embedding (VE) is incorporated into the FPX code to enhance the prediction 
capability of parallel blade- vortex interactions 11 . 

A grid generation package GRGN3 is used to generate 0-H mesh around rotor blade. 
The blade surface is defined by an input file. The far-field boundary of the mesh is set at a 
fixed number of chords from blade surface. The mesh points are generated between blade 
surface and far-filed boundary. This grid generator is combined into the flow solver FPX 
recently 12 . 


4. Implementation 

The coupling of the CDISC design method with the rotor FPX code is the main scope 
of the present work. Two methods can be used for coupling the CDISC module with 
the FPX flow analysis code: direct coupling where the CDISC module is treated as a 
subroutine called by the flow solver; and indirect coupling where a job control language is 
used. The indirect coupling offers advantages in code maintenance and portability, and the 
design module does not add to the central memory requirement. In the present work, an 
indirect coupling is used where a UNIX script file is used to control entire design procedure 
automatically. User interference during the design processes is not needed. 

Information is passed between the CDISC module and the FPX flow solver using blade 
section geometry file and surface pressure coefficient file. The FPX flow solver with a built- 
in grid generator 12 is used in this coupling, so that the grid can be re-generated within the 
flow solver each time when a new blade section geometry file is obtained from the CDISC 
module. Some modifications on input and output statements to the CDISC code and the 
FPX code are made as required by the coupling, and thus any type of auxiliary pre- or 
post-processor codes is not needed. The UNIX script file starts the FPX flow solver to 
produce an output file containing blade section surface pressure coefficients. This file along 
with the initial blade section geometry file is then fed into the CDISC design module. The 
CDISC design module then produces a new geometry file. This new geometry file is used 
as input for the FPX flow solver. The built-in grid generator generates anew grid and then 
the flow solver produces a new surface pressure coefficient file. This procedure continues 
for desired number of cycles as specified in the UNIX script file. Typical cases take 15-30 
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cycles to get a convergent design. 


5. Results 

As a test of capability of the new package using the CDISC and the FPX codes in 
rotor section design, several sample design cases are made for a rotor blade. The planform 
of the blade is shown in Figure 1 along with seven design stations as indicated. The 
blade has a zero-twist. The first three cases are basic DISC design where target pressure 
distributions are pre-defined without any constraints. The next two cases are CDISC 
design where target pressure distributions are automatically generated. The following 
sub-sections present these results. 

5.1 Basic DISC Results 

In all three basic DISC design cases, initial blade sections are NACA0012 airfoil sec- 
tions and target pressure distributions are obtained from a known target geometry. The 
new design package is then used to “reproduce the original target configuration. The 
tip Mach number is set at 0.6288 for the first and second cases, 0.95 for the third case. 

In the first case, the “target” geometry is NACA0008 airfoil section. Thirty design 
cycles are executed to get a convergent (almost convergent at tip-station) design. This case 
requires around 26 minutes on the Cray-C90 computer with a single processor. The results 
of this case are presented in Figures 2a-2e. Figures 2a and 2b are design history in terms 
of blade section geometry and surface pressure distributions, respectively. Figure 2c gives 
design history of spanwise lifting coefficients, moment coefficients and drag coefficients. 
A good approximation to the target pressure distributions is achieved after only a few 
(ten) design cycles (target pressure is given next). Figures 2d gives a comparison of design 
geometry and target geometry along with the initial geometry, while Figure 2e gives a 
comparison for surface pressure distributions for four design stations as indicated (although 
a total of seven stations are designed). The correlation between the final design and target 
pressure distributions at each station is very good. The new airfoil sections generated 
compare well with the original “target” geometry. 

The second case is very similar to the first case. In this case however, the “target” 
geometry is an “after- design 1 NACA0008 airfoil, where a CDISC procedure is first applied 
to NACA0008 to produce a modified airfoil for reducing drag. The target pressure dis- 
tribution is the final design pressure during this design process. The results of this case 
are presented in Figures 3a-3e. Figures 3a-3c are design history. Figures 3d and 3e are 
comparison of design and target. A good agreement between design and target is obtained 
again. 

In the third case, the “target” geometry and target pressure are provided from a 
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CDISC design case applied to MDHC HH02 and NACA64A006 airfoil sections. It is ex- 
pected that a large number of design cycles is needed for this case to get a convergent 
design due to a substantial difference between the initial geometry and the “target” ge- 
ometry, particularly at tip station. As a consequence, 300 design cycles are executed for 
this case to examine the nature of convergence. The results of this case are presented in 
Figures 4a-4g. Figures 4a and 4b are design history in terms of airfoil section geometry 
and surface pressure distributions, respectively. Figure 4c gives design history of spanwise 
lifting coefficients, moment coefficients and drag coefficients. It is seen that the solution 
is convergent between design cycles 150 and 200 for Stations 1, 3 and 5, and almost con- 
vergent at design cycle 300 at Station 7 (the tip station). Figures 4d gives a comparison 
of final design geometry and target geometry along with the initial geometry, while Figure 
4e gives a comparison for surface pressure distributions for the four design stations. The 
correlation between the final design and target pressure distributions at Stations 5 and 7 
is very good. However, some discrepancies exist for Stations 1 and 3. Same is true for the 
comparison of geometry. Figures 4f and 4g then made another comparison of target with 
design where design geometry and pressure are obtained at design cycle 50 for Stations 1, 
3 and 5. It is now noted that the agreement between design and target is better than that 
of Figures 4d and 4e for Stations 1 and 3, particularly for Station 1. 

5.2 CDISC Results 

With the Constrained DISC approach, the emphasis is shifted, to meeting the con- 
straints where target pressure is generated automatically within the design module. In the 
next sample cases, two constraints are imposed: one to fix the original spanload distri- 
bution, and other one to reduce drag using the Modified Uniform Distribution (MUD) of 
chordwise loading constraint option. A MUD value of 0.2 is used according to Campbell 5 
to obtain a best result in reducing drag. The tip Mach number is set at 0.95. In both 
cases, 15 design cycles are executed. 

Figures 5a-5c are the results of the first CDISC sample design case, where initial blade 
sections are NACA0012 sections. Figures 5a and 5b give the history of design for blade 
section geometry and surface pressure distributions, respectively. Figure 5c presents a 
spanwise lifting coefficient, moment coefficient, drag coefficient and maximum local Mach 
number distributions. It is seen that the design reduces the negative lifting force, pitching 
moment and drag force. By examining the drag coefficients, it is found that the drag is 
reduced by 12.5% at Station 3 and 7.5% at Station 5 through the CDISC design. The 
maximum local Mach numbers are unchanged during the design process. 

Figures 6a-6c are the results of the second CDISC sample design case, where initial 
blade sections are MDHC HH02 airfoil for Stations 1 through 5, and NACA64A006 airfoil 
for Stations 6 and 7. Figures 6a and 6b give the history of design for blade section geom- 
etry and surface pressure distributions, respectively. Figure 6c presents a spanwise lifting 
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coefficient, moment coefficient, drag coefficient and shock Mach number distributions. By 
examining the drag coefficients, it is found that the drag is reduced by 12.7% at Station 1, 
9.6% at Station 3, 25.6% at Station 5 and 12.5% at Station 7 through the redistribution 
of the loading and smoothing. Shock is also eliminated through the design as indicated, 
in the figure, by zero-shock Mach number along the span. 

6. Conclusions 

A rotor section aerodynamic design package is developed by coupling the CDISC 
design method with the FPX rotor CFD code. The coupling between the CDISC design 
and the FPX flow analysis is fully automated using a UNIX script file to control job flow. 
Five sample design cases are presented for both basic DISC and Constrained DISC design 
for a zero-twist blade in hovering motion. The results indicate that the CDISC method is 
a robust and efficient design method when used in rotor section design coupled with the 
FPX code. 

These sample cases demonstrate that the target pressure can be achieved even when 
the initial geometry is not “close” to the “target” geometry. A constrained DISC design 
case demonstrates a shock-free design using MUD method, where drag is significantly 
reduced. Actual computational effort depends on the complexity of the target pressure 
distributions and the difference between the initial and “target” geometry. A typical 
design case with 30 design cycles for designing 7 stations requires around 26 minutes CPU 
time on the Cray-C90 computer with a single processor. Presently, all sample design cases 
are made to a zero-twist blade without in-flow angles for either symmetric or cambered 
airfoil section design. The further work is needed in extending the current design package 
into high-lifting flows with large twist and in-flow ang les 
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Figure 1. The rotor blade planform. 
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Figure 2a. Design history of blade section geometry 
for basic DISC design case 1, Mu p == 0.6288. 
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Figure 2b. Design history of surface pressure distributions 
for basic DISC design case 1, M t i p = 0.6288. 
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coefficients for basic DISC design case 1, M tip = 0.6288. 
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Figure 2e. Comparison of surface pressure distributions 
for basic DISC design case 1, M ttp = 0.6288. 
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Figure 3a. Design history of blade section geometry 
for basic DISC design case 2, M tiv = 0.6288. 








Figure 3b. Design history of surface pressure distributions 
for basic DISC design case 2, M txp — 0.6288. 
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Figure 4a. Design history of blade section geometry 
for basic DISC design case 3, M tip = 0.95. 
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Figure 4b. Design history of surface pressure distributions 
for basic DISC design case 3, M ttp = 0.95. 
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Figure 4d. Comparison of blade section geometry for basic DISC 
design case 3, M t i p — 0.95, where final design is used. 
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Figure 4e. Comparison of surface pressure distributions for basic DISC 
design case 3, M ttp = 0.95, where final design is used. 
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Figure 4f. Comparison of blade section geometry for basic DISC 
design case 3, Mttp = 0.95, where the design at cycle 50 is used. 
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Figure 5b. Design history of surface pressure distributions 
for CDISC design case 1, M t i p = 0.95. 
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r/R (in %) r/R (in %) 

Figure 5c. Design history of spanwise lifting ( C L ), moment (Cm), drag (Co) 
coefficients and local maximum Mach No. for CDISC design case 1, M tiv = 0.95. 
















Figure 6b. Design history of surface pressure distributions 
for CDISC design case 2, M tip = 0.95. 
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Figure 6c. Design history of spanwise L 
coefficients and shock Mach No. for 
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